Molecular dynamics study of contact mechanics: contact area and interfacial 

separation from small to full contact 



o 
o 

(N 
-t— > 

o 
O 



o 

s 

I 

c 

o 
o 



> 

l> 
(N 
O 

d 
o 



X 



Oh 

to 

O 



C. Yang and B.N.J. Persson 
IFF, FZ-Julich, 52425 Jiilich, Germany 

We report a molecular dynamics study of the contact between a rigid solid with a randomly 
rough surface and an elastic block with a flat surface. We study the contact area and the interfacial 
separation from small contact (low load) to full contact (high load). For small load the contact area 
varies linearly with the load and the interfacial separation depends logarithmically on the load. For 
high load the contact area approaches to the nominal contact area (i.e., complete contact), and the 
interfacial separation approaches to zero. The present results may be very important for soft solids, 
e.g., rubber, or for very smooth surfaces, where complete contact can be reached at moderate high 
loads without plastic deformation of the solids. 



It is a very difficult to prepare surfaces which are really 
flat. Even on most polished surfaces, hills and valleys are 
present which are large compared with the atomic-size. 
Usually, if two solids are placed in contact, the upper 
surface will be supported on the summits of the irreg- 
ularities, and large surface areas will be separated by 
distances which are great compared with the molecular 
range of action[l|, 0, S S|- The separation it(x) between 
the surfaces will vary in a nearly random way with the 
lateral coordinates x = (x, y) in the apparent contact 
area. When the applied squeezing pressure increases, the 
contact area A will increase and the average surface sep- 
aration u = (u(x)) will decrease, but in most situations 
it is not possible to squeeze the solids into perfect con- 
tact corresponding to u = 0. Understanding the area of 
real contact, and the interfacial separation between two 
solids is essential to friction, adhesion, sealing and many 
other important applications 0, Q. 

Most studies of contact mechanics have been focused 
on small load where the contact area depends linearly 
on the load@, H, [1, 03, 11 1. However, for soft solids, 
such as rubber or gelatin, or for smooth surfaces, nearly 
full contact may occur at the interface, so it is of great 
interest to study how the contact area, the interfacial 
surface separation and stress distribution vary with load 
from small load (where the contact area varies linearly 
with the load) , to high load [where the contact is (nearly) 
complete] . Here we will present such a study using molec- 
ular dynamics(MD), and we will compare the numerical 
results with the prediction of the analytical contact me- 
chanics theory of Persson 12,13, [13] • 
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Consider randomly rough surfaces with roughness 
wavelength components in some finite range Ai < A < Ao, 
where Ao is similar to (but smaller than) the lateral size 
of the nominal contact area. In order to accurately de- 
scribe the contact mechanics between elastic blocks, it is 
necessary to consider solid block which extends (at least) 
a distance ~ Ao in the direction normal to the nominal 
contact area. This leads to an enormous number of atoms 
or dynamical variables even for a small systems. In order 
to avoid this trouble we have developed a multiscale MD 



FIG. 1: The pressure distribution for £ = 4 for three differ- 
ent nominal pressure, (a) p/E* w 0.013, (b) p/E* « 0.113, 
(c) p/E* ~ 0.218. The pressure probability distribution be- 
comes broader with increasing the squeezing pressure. Prop- 
erly choosing the correction factor (see Fig. [2} makes the 
numerical results in good agreement with Persson's contact 
mechanics theory. 



approach [13]. The atoms at the interface between the 
block and substrate interact with the repulsive potential 
U(r) = e (r /r) 12 , where r is the distance between a pair 
of atoms, r = 3.28 A and e = 74.4 meV. In the MD- 
model calculations there is no unique way to define the 
separation u between the solid walls. Here we have used 
the same definition as in Ref. jlQI ] u = d — d c , where d 
is the difference between the plane through the center of 




o 

a 

o 
u 

o 
U 



0.2 0.4 0.6 0.8 
A/Ao 

FIG. 2: (a) Correction factor r as a function of the contact 
area ratio A/Ao- The points are simulation results, which 
have been fitted by the function f(x) = a + bx 2 + cx 3 + dx 4 
under the condition d — 1 — a — b — c (solid line). The corre- 
sponding coefficients a, b, c are 0.51, 2.5, —3.3 respectively, 
(b) The contact area ratio A/Ao as a function of squeezing 
pressure normalized by effective elastic modulus. 



the atoms of the top layer of substrate atoms and bottom 
layer of block atoms. d c is the critical atom-atom sepa- 
ration used to define contact on the atomic scale. Thus, 
u = corresponds to the separation d c = 4.36A between 
planes through the center of the interfacial atoms of the 
block and the substrate. 

The system has lateral dimension L x = N x a and 
L y = N y a, where a is the lattice space of the block. 
Periodic boundary condition is used in xy plane. For 



N„ 



400, while the lattice space of 



the block N x 

the substrate b » a/</>, where <f) = (1 4- \fh)/2 is the 
golden mean, in order to avoid the formation of com- 
mensurate structures at the interface. The mass of the 
block atoms is 197 a.m.u. and a — 2.6 A, reproducing 
the atomic mass and density of gold. The elastic mod- 
ulus and Poisson ratio of the block is E = 77.2 GPa 
and v = 0.42, respectively. The substrate surface has 
self-affine fractal surface roughness (Tol Il8j |. For a self- 
affine fractal surface the power spectrum has power-law 
behavior C(q) ~ q~ 2 ( H + 1 ) ; where the Hurst exponent 
H is related to fractal dimension Df of the surface via 
H = 3— Df. For real surfaces this relation holds only for a 
finite wave vector region qo < q < q\ . Note that in many 
cases there is roll-off wavevector qo below which C(q) is 
approximately constant. The randomly rough substrates 
we use have been generated as described in Ref . [3] , and 
have root-mean-square roughness h rms — 10A, fractal di- 
mension is Df = 2.2, and roll-off wavevector qo = 3ql, 
where q^ — 2tt/L. We define the magnification £ = q/qo- 

With molecular dynamics simulations we can calculate 
the interfacial stress distribution. In order to obtain the 
contact area we follow the procedure outlined in Ref. 
and fit the numerical results to the theoretically predicted 



stress distribution 
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where G(p, C) depends on the nominal squeezing pressure 
p and the magnification £ (but which is independent of a, 
see below), and which we choose to get the best fit with 
the numerical data. In Fig. [1] we have shown the good 
agreement between the numerical pressure distribution 
and the analytical theory (for £ = 4) under three different 
nominal pressure. Once G is known we can calculate the 
relative contact area using 0] 
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In Fig. HJb) we show the relative contact area A/A$ as a 
function of normalized pressure p/E* from small to full 
contact. 

The fitted G(p, Q can now be compared with theory. 
Thus, the theory of Persson predicts G = G where 
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Fig. [2ja) shows the ratio r — G/G which refers to cor- 
rection factor. Note that r increases from sa 0.51 to 1 as 
the squeezing pressure p increases from zero to infinite 
(i.e., the normalized contact area A/Aq increases from 
zero to 1). Since the contact area for small load is pro- 
portional to ~ 1/%/G, it follows that the theory for small 
load predicts a contact area about ~ 30% smaller than 
that deduced from the MD simulation. This is slightly 
larger than what has been found in earlier numerical sim- 
ulations. Thus, the finite element calculations of Hyun 
and Robbins[26j and the Green's function molecular dy- 
namics study of Campana and Miiser 11] gives r « 0.64, 
corresponding to a contact area about ~ 20% larger than 
that predicted by the Persson theory. Similarly, the study 
of H6nig[27| gives r i=s 0.56 for small load. However, 
none of the computer simulations can be considered as 
perfectly converged, so the difference between theory and 
fully converged numerical simulation may be smaller than 
that indicated by the numbers given above. 

Recently Persson theoretically derived the relation be- 
tween the average interfacial separation u and the applied 
normal squeezing pressure p [141 ] . For nonadhesive inter- 
action and small applied pressure, p ~ e -«/«o j m a good 
agreement with recent experimental observations [22| . 
Here we numerically calculate the average interfacial sep- 
aration with different squeezing pressure with molecular 
dynamics. In Fig. [3]we show the natural logarithm of the 
normalized average pressure p/E*, as a function of the 
normalized interfacial separation u/h rms . We show result 
for the magnification £ = 4 (open circles) and £ = 216 
(solid squares). Since the atoms interact with a long- 
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FIG. 3: An elastic block squeezed against a rigid rough sub- 
strate. The natural logarithmically normalized average pres- 
sure log(p/E*), as a function of normalized interfacial sepa- 
ration u/h ma on different magnifications £ = 4 and £ = 2f6. 





FIG. 5: Probability distribution of interfacial separation un- 
der different pressure (a) p/E* « 0.002, (b) p/E* w 0.013, (c) 
p/-E* « 0.II3, at low and high magnifications respectively. 
Note that the distributions of interfacial separations observed 
at low and high magnifications are similar for u > 5 A. This 
result is expected since mainly the long- wavelength, large am- 
plitude roughness will determine the separation between the 
surfaces when the separation is large. 



FIG. 4: The relation between the natural logarithm of the 
squeezing pressure p (normalized by E*) and the interfacial 
separation u (normalized by the root-mean-square roughness 
amplitude farms) for an elastic solid squeezed against a rigid 
surface. The theory curve (solid line) has been calculated 
using the theory presented in Ref. [l4| with 7 = 0.42. 

squeeze the surfaces closer to each other than what cor- 
responds to u = 0. This explains why simulation data 
points occur also for u < 0. 

In Fig. @]we compare the MD results from Fig. [3] (open 
circles) with the theory presented in Ref. [3] using the 
same surface roughness power spectra (and other param- 
eters) as in the MD-calculation. The theory is in good 
agreement with the numerical data for 0.2 < u/h rms < 2. 
For u/h lma < 0.2 the two curves differ because of the 
reason discussed above, i.e., the "soft" potential used in 
the MD simulation allows the block and substrate atoms 
to approach each other beyond u(x) — 0, while in the an- 
alytical theory the potential is infinite for w(x) < and 
zero for u(x) > 0. The difference between the theory and 
the MD results for u/h Yms > 2 is due to a finite size ef- 



fect. That is, since the MD calculations use a very small 
system, the highest asperities are only ~ 3/i rms above the 
average plane (see the height distribution in Ref. [29|). 
and for large u very few contact spots will occur, and in 
particular for u > 3/i rms no contact occurs and p must 
vanish. In the analytical theory, the system size is as- 
sumed to be infinite. Even for a Gaussian distribution of 
asperity height, there will always be infinitely many in- 
finitely high asperities. Contact will occur at arbitrarily 
large separation u, and the asymptotic relation u ~ logp 
will hold for arbitrarily large u at small squeezing pres- 
sures p. 

Let us now discuss the probability distribution of in- 
terfacial separation, defined by 

P u = (6(u u(x))) (4) 

where (...) is ensemble average. This function is shown in 
Fig. [5] for three different loads and two different magnifi- 
cations. Note that the distributions of interfacial separa- 
tions observed at low and high magnifications are similar 
for u > 5 A. This result is expected since mainly the long- 
wavelength, large amplitude roughness will determine the 



separation between the surfaces when the separation is 
large. The quantity P u has many important applications. 
For example, for lubricated contact at low sliding veloc- 
ity, one may estimate the contribution from shearing the 
liquid film to the (nominal) frictional stress using 

f°° Pu 
a sa nv / du — - 

Ju c u 

where r/ is the viscosity and v the sliding velocity, and 
where u c is a cut-off separation of order nanometer (con- 
tinuum fluid dynamics is not valid for liquid films thinner 
than a few nanometers). Another important application 
is for estimating the fluid leaking through sealing[29]. 

To summarize, we have performed a Molecular Dy- 
namics (MD) study of the contact between an elastic 
block with a flat surface and a rigid substrate with a 
randomly rough surface. The interfacial pressure distri- 
bution agrees well with the analytical theory of Persson. 
We have also calculated the area of real contact and the 
interfacial separation as a function of load from small 
to full contact, and compared the results with the an- 
alytical theory. For not too large or too small squeez- 
ing pressures, the MD results show that the interfacial 
separation u depends logarithmically on the squeezing 
pressure u ~ logp, in a good agreement with the analyt- 
ical theory and with experimental observation. Finally, 
we have studied the probability distribution of interfacial 
separation for different pressure and different magnifica- 
tions. The present results may be of great importance 
for soft solids, e.g. rubber-like material, or very smooth 
surfaces. 
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